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We investigate the magnetocaloric properties of the two-dimensional frustrated J1-J2 model on 
a square lattice. This model describes well the magnetic behavior of two classes of quasi-two- 
dimensional S = 1/2 vanadates, namely the Li2VOX04 (X = Si, Ge) and y4y4'VO(P04)2 {A, A' = 
Pb, Zn, Sr, Ba) compounds. The magnetocaloric effect (MCE) consists in the adiabatic temperature 
change upon changing the external magnetic field. In frustrated systems, the MCE can be enhanced 
close to the saturation field because of massive degeneracies among low lying excitations. We discuss 
results for the MCE in the two distinct antiferromagnetic regimes of the phase diagram. Numerical 
finite temperature Lanczos as well as analytical methods based on the spin wave expansion are 
employed and results are compared. We give explicit values for the saturation fields of the vanadium 
compounds. We predict that at subcritical fields there is first a (positive) maximum followed by 
sign change of the MCE, characteristic of all magnetically ordered phases. 

PACS numbers: 75.10.J, 75.40.C 
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I. INTRODUCTION 

Two dimensional (2D) magnetic systems are favorite 
models to study the influence of quantum fluctuations on 
magnetic order. Depending on the model they may both 
prohibit an ordered ground state or select a specific order 
among classically degenerate states. These phenomena 
have been studied in great detail for geometrically frus- 
trated systems like trigonal, Kagome and checkerboard 
lattice^. However they are also present in magnets where 
the frustration is not the result of lattice geometry but of 
competition between different (for example nearest- and 
next-nearest neighbor) magnetic bonds. A prime exam- 
ple is the frustrated J1-J2 model on a square lattice. Its 
ground state and thermodynamic properties in zero field 
have been well studied (see Refs. [HHH and references 
cited therein). 

Classically the model predicts three magnetic phases 
depending on the frustration ratio J2/J1: The ferromag- 
net (FM), (tt, tt) Neel antiferromagnet (NAF) and (tt, 0) 
coUinear antiferromagnet (CAF). However it is known 
that close to the classical CAF/NAF and CAF/FM 
boundary quantum fluctuations destroy magnetic order 
and presumably stabilize nonmagnetic order parameters. 

The discovery of two classes of layered vanadium ox- 
ides Li2VOX04 (X = 31,06)^ and y4y4'VO(P04)2 
{A, A' = Pb, Zn, Sr, Ba.)^ which are weU described by 
this model has further raised interest in the J1-J2 model. 
One advantage of the new vanadium compounds is a com- 
paratively low energy scale for the exchange constants of 
order 10 K. Therefore high field experiments might be a 
promising way to learn more about their physical prop- 
erties, indeed the saturation field for these compounds 
where the fully polarized state is achieved seems within 
experimental reach. 

Therefore in this work we study exhaustively the high- 
field magnetic and especially the magnetocaloric effects 



(MCE) in the J1-J2 model. We use a variety of ana- 
lytical and numerical techniques to investigate the de- 
pendence of magnetization, susceptibility, entropy spe- 
cific heat and adiabatic cooling rate on magnetic field, 
temperature and frustration ratio. The variation of the 
saturation field with the frustration ratio is calculated 
and predictions for the abovementioned compounds are 
made. We show that the low temperature specific heat is 
strongly enhanced around the classical phase boundaries 
where large quantum fluctuations occur. 



Our special focus is on the magnetocaloric effect. We 
will show that the cooling rate, normalized to its para- 
magnetic value, is strongly enhanced above the satura- 
tion field and depends on the frustration angle. We also 
predict that for subcritical fields the cooling rate is first 
positive with a maximum at moderate fields and a sub- 
sequently changes sign at a larger subcritical field. This 
behavior is common to all AF phases of the model and 
can be understood quantitatively from calculations of 
contours of constant entropy in the (/i, T) plane. The 
dependence of corresponding characteristic fields on the 
frustration ratio are also calculated. 



In Sec. |TT] we give a brief description of the magne- 
tocaloric effect in magnets. In Sec. Illll the basic proper- 
ties and phase diagram of the J1-J2 model are introduced. 
In Sec. lIVI we discuss extensively results of the finite tem- 
perature Lanczos method (FTLM) for finite 2D J1-J2 
clusters. In Sec. |V] we use analytical methods within 
mean field or spin wave approximation as an alternative 
way to study the magnetocaloric properties. In Sec lVIl 
we discuss and compare the results obtained by the var- 
ious methods. Finally Sec. IVIll gives the summary and 
conclusion. 
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II. THE MAGNETOCALORIC EFFECT IN 
MAGNETICALLY ORDERED COMPOUNDS 

When a crystal containing magnetic ions is placed in a 
magnetic field the adiabatic or isentropic change of this 
external parameter causes a temperature change in the 
sample. This is called the magnetocaloric effect (MCE) 
which was first discovered by Warburg^. It is nowa- 
days interesting in several aspects. Firstly suitable com- 
pounds, like paramagnetic salts where demagnetization 
leads to cooling may be used technicallyi^iii. Secondly 
at high (pulsed) fields the magnetocaloric anomalies at 
a magnetic phase transition may be used to map out 
the H-T phase diagrams which are not accessible oth- 
erwise. Finally it has recently gained special attention 
in frustrated magnets. There the behavior around the 
saturation field may be described by the condensation 
of a macroscopic number of local magnons^^'^^'^"' which 
leads to a giant enhancement of the magnetocaloric cool- 
ing rate. The latter is defined as the rate of change of 
temperature T with magnetic field H at fixed entropy S. 
Using a Maxwell relation one can write this as 



dT 
'dH 



Cv V dT 



\dT) H 



(1) 



H 



where Cy is the heat capacity and M the magnetiza- 
tion of the sample. The integrated adiabatic tempera- 
ture change along an isentropic line with S{T, H) — const 
which is caused by the variation of magnetic field is then 
given by 

AT,d{Ho,H) = T-To= / r,,,,{H' ,T')dH' (2) 

Jh„ 

Here Tq and Hq are the starting values of temperature 
and field respectively. We take the adiabatic cooling rate 
for free paramagnetic ions as a reference quantity. As 
shown later it is simply given by = (T/H). The 
dimensionless magnetocaloric enhancement factor due to 
interaction effects is then defined by the ratio 



f mc = ^mc/^mc ^ {H/T)r„ 



(3) 



In the present work we study the MCE on a Heisenberg 
square lattice J1-J2 model which incorporates a frustra- 
tion of nearest- and next-nearest neighbor exchange in- 
teractions. It is a suitable model to analyze the mag- 
netic properties of two classes of quasi-two-dimensional 
vanadates, namely the Li2VOX04 {X = Si, Ge)ii^i^ and 
^A'VO(P04)2 {A, A' = Pb,Zn,Sr,Ba) compoundsli^. 
We investigate several aspects of magnetocaloric prop- 
erties. We show that indeed it may be used to iden- 
tify the saturation field and the associated Hc[T) phase 
boundary between the fully polarized and the AF ordered 
states. We also discuss whether it may be used as a diag- 
nostic for the appropriate frustration angle (or J1I J\ ra- 
tio). Finally we will study whether the frustration effect 



lead to a visible signature in the anomalies of the magne- 
tocaloric cooling rates, especially close to the saturation 
fields. We will employ both numerical FTLM methods as 
well as approximate analytical methods based on mean 
field and spin wave approximations. 



III. THE J1-J2 HEISENBERG MODEL AND 
EXAMPLES 

We first give a brief characterization of the 2D square 
lattice J\-Ji model in an external field which is defined 
by the Hamiltonian 



n 



Ji + J2 ^ S,.S, -hY^Sl (4) 



(y>i 



{'y>2 



with the convention that each bond is counted only once. 
Here J\ is the nearest-neighbor exchange coupling along 
the edges and J2 the next-nearest neighbor exchange cou- 
pling along the diagonals of each square. Furthermore 
h = g^iBH where H is the applied magnetic field. Here 
g is the gyromagnetic ratio and /xb the Bohr magneton. 
The zero-field phase diagram may best be characterized 
by introducing the equivalent parameter set 



Jc = [J^ + J^)i; 4> = tan-i(J2/Ji) 



(5) 



The 'frustration angle' is a convenient quantity to 
characterize the amount of exchange frustration in the 
model and Jc gives the energy scale at which thermody- 
namic anomalies in specific heat susceptibility etc. are to 
be expected. For spin-1/2, as function of (j) three main 
phases (FM,NAF,CAF) appear already on the classical 
level (Fig. However for <j) ~ O.IStt (J2/ Ji ~ 0.5) 

and (j) ~ 0.857r ( J2/ Ji 0.5) where CAE meets NAF 

and FM respectively a large degeneracy of the classical 
ground state appears and quantum fluctuations lead to 
non-magnetic "hidden order" phases shown as the shaded 
sectors in Fig. [1] These phases have been extensively dis- 
cussed iniii2ii^ and references cited therein. 

For the compounds mentioned above generally 
Jc ~ 10 K. In addition from the high temperature expan- 
sion of the susceptibility the Curie- Weiss temperature is 
obtained as Qcw = {Ji + J2)/fcB. Both quantities can 
be obtained from experiment and the appropriate pair of 
exchange constants (Ji, J2) then has to lie on the inter- 
sect of a circle ( Jc) and a straight fine {&cw) as shown 
in Fig. [TJ Obviously there are always two solutions lying 
in the NAF (</>-) and CAE {(t>+) sector. This observation 
is unchanged by a more detailed analysis of susceptibility 
and specific heat^ii^. 

Various other methods have been proposed to resolve 
the ambiguity of frustration angles such as measurement 
of the spin structure factor—, nonlinear susceptibilityi^ 
and the saturation field of the magnetizatior>ii. Only 
the former has been tried sofar for Li2VO(Si04) and 
Pb2VO(P04)5r^. In both cases the ground state clearly 
has CAE order. For this reason we have assigned the 
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J2/J1= -00 



FIG. 1: Phases of the spin-1/2 2D square lattice J1-J2 modeL 
The CAP and NAP order is indicated by arrows, the associ- 
ated wave vectors are Q — (0,1) or (1,0) and (1,1) (in units 
of 7r/a) respectively. The boundary between PM and NAP 
phase is the line Ji = 0, J2 < 0. Values of J2/J1 in parenthe- 
ses indicate where zero point fluctuations destroy the CAP 
order parametei^. Dashed lines correspond to the experimen- 
tal 0CW = {Ji + J2)/k-B and refer to the known J1-J2 com- 
poundsii^. Two solutions 4>+ (CAP) and 4>- (NAP) are com- 
patible with the thermodynamic properties. Here we choose 
the former since they are confirmed by neutron diffraction for 
the Li- and Pb- compounds 



other known family members of J1-J2 vanadates to the 
same sector in Fig.[T]although a confirmation for this con- 
jecture is still lacking. We believe that high field investi- 
gations are a further promising method to shed light on 
these compounds , especially because the comparatively 
low Jc (~ 10 K) will lead to saturation fields relatively 
easy to access. 

Therefore in this work we study the J1-J2 model in 
an external field as given by Eq. Thereby we fo- 

cus on the theory of the saturation field and the magne- 
tocaloric anomalies both around the saturation field and 
for smaller fields within the ordered phase. We will use 
both numerical analysis of finite clusters based on the 
FTLM method as well as analytical methods based on 
mean field or linear spin wave approximations for com- 
parison. In the latter we focus on the three main phases 
with magnetic order. The analytical treatment of the 
magnetocaloric effect in the hidden order phases war- 
rants a separate treatment which takes into account the 
proper non-magnetic order parameter. 



IV. EXACT DIAGONALIZATION FOR FINITE 
CLUSTERS AT FINITE TEMPERATURES AND 
FINITE FIELD 

We have performed numerical exact-diagonalization 
calculations for three different clusters; squares with 



16 and 20 sites and a 24-site rectangle. All of these 
tile the lattice in such a way as to be compatible with 
both the (tt, tt) NAF and(7r, 0) CAF states, once periodic 
boundary conditions are imposed. Our main focus is on 
the finite-temperature, finite-field properties of the J1-J2 
model. We therefore use the finite-temperature Lanc- 
zos method (FTLM) to evaluate the partition function 
of the model, together with thermodynamic averages of 
the form 

{A{T,H)) = lTr(Ae-«/('=BT))^ (g) 

Z = Tre-^/C^-^), (7) 

Here A is an operator, Ti is the Hamiltonian of the J1-J2 
model including the Zeeman term (Eq. ^) and 2(T, H) 
its (field-dependent) partition function. For each frustra- 
tion angle 0, we perform between 100 and 500 Lanczos 
iterations with different starting vectors in each symme- 
try sector of the Hilbert space. We keep between 1 and 
100 eigenvalues and eigenvectors of the tridiagonal Lanc- 
zos matrix per iteration in order to evaluate the thermo- 
dynamic traces discussed in the following. Details of the 
method can be found in Ref. il9i . 



A. Level crossings, spinwave instabilities and 
saturation fields 

Before addressing the finite-temperature results, let us 
examine certain general features of the model at zero 
temperature. Applying a magnetic field H leads to 
a Zeeman splitting of the energy levels, and therefore 
level crossings occur when increasing the field. These 
level crossings correspond to jumps in the magnetiza- 
tion at zero temperature, until the fully polarized state 
is reached at a certain critical value of the magnetic field. 
Figure [2] illustrates this behavior for two different values 
of the frustration angle </>: 

As an example for positive (antiferromagnetic) Ji, we 
show the field dependence of the energy levels for the 
pure Neel antiferromagnet (J2 = 0) on the left side of 
the figure. In the whole "right half" of the phase dia- 
gram ( Ji > 0, J2 arbitrary), the field dependence of the 
energy levels is qualitatively similar. For the 24-site clus- 
ter considered, 12 spin flips with AS^ = 1, indicated by 
the small arrows, occur at the points where the magnetic 
field is given by 

ffMBi?flip(5.) = ^ (So(5.) - £^0(5. - 1)) . (8) 

Here, Eo{Sz) denotes the ground-state energy for the 
subspace with constant Sz at zero field where N is the 
cluster size. The saturation field Hgat is reached when 
Sz = N/2. Because the fully polarized state is an eigen- 
state of the Hamiltonian, the numerical values for -ffgat 
from the equation above are exactly identical to what 
one finds within linear-spinwave theory for the infinite 
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FIG. 2: Energy levels as a function of the applied magnetic field for two different frustration angles </!> = (Neel antiferromagnet, 
Ji > 0, left) and ^/tt = 0.84 (coUinear antiferromagnet, Ji < 0, right). For each sector of the Hilbert space with total 
Sz = const., the field dependence of the respective ground state is plotted. The arrows point to the energy/field values where 
a jump in the ground-state magnetization of the full system occurs. 



system, 



5MB rrLSW 
T -"sat 



zS 



cos <j)[l — - (cos Qx + cos Qy) 



+ sin 1 — cos Qx cos ( 



(9) 



with z = A, S = 1/2, and Q = (tt, tt) or one of (tt, 0), 
(0, tt) is the antiferromagnetic ordering vector. 

In contrast, for ferromagnetic Ji < (but outside 
the ferromagnetic regime in the phase diagram) , we see a 
quahtatively different behavior of the energy levels. The 
right-hand-side of Fig. [2] shows an example for the frus- 
tration angle ^/tt = 0.84. Instead of 12 level crossings, 
there are only six, each corresponding to a jump \Sz = 2 
occurring at fields 



1 

iVfc 



{Eo{S,)^E^{S,~k)), fc-2. (10) 



The saturation field is now given by an instability crite- 
rion of the fully polarized state towards a two-magnon 
excitation, 

(11) 

For Ji < 0, this field is larger than the field of the one- 
magnon instability given by the above equation with k = 
1 and therefore determines the predominant instability 
when lowering the field in the fully polarized state. 

A necessary condition for a AS'^ = 1 level crossing 
to occur is that the lower bound Eq{Sz) of the energy 
spectrum at zero field for a fixed value Sz is a convex 
function of Sz, i.e., the condition 



Eo{Sz + l) <l{Eo{S, 



)+Eo{Sz + 2)) (12) 



must be fulfilled at H = 0. At the special point Ji = 0, 
J2 > (0 = 7r/4), the J1-J2 lattice decouples into two 



independent Neel sublattices, and equality holds above. 
For the finite size clusters which we consider, enlarg- 
ing (j) further (i. e. making Ji ferromagnetic) stabilizes 
two-magnon bound states, and Eq. no longer holds. 
Level crossings are characterized by AS*^ — 2, and the 
saturation field iJc2 is given by Eq. ([TT]) with k = 2. 
This is exactly what would be expected where a spin ne- 
matic state is selected by quantum fiuctuations in applied 
magnetic field^. 

However these are finite size results, and must be ap- 
proached with a little caution. The critical fields asso- 
ciated with one- and two-magnon excitations show quite 
different finite size scaling properties as a function of 0, 
as illustrated in Fig. [31 From the three different clus- 
ter sizes studied here, the following observations can be 
made : Firstly, AHc is a non-monotonic function of the 
frustration angle; it has a minimum AHc = at the 
crossover between the NAF and CAE phases for Ji — 2J2 
{(j>/n fa 0.15). Secondly, AHc changes sign at Ji = in 
the CAE phase {4>/Tr = 1/2) in favor of a ASz = 2 insta- 
bifity as above. Thirdly, for -1/2 < ^/tt < 0.8, \AHc\ 
is a monotonically decreasing function of 1/iV and seems 
to extrapolate to zero for N ^ 00. This means that the 
one-magnon instability (and conventional canted AF or- 
der) is restored in the thermodynamic limit for FM Ji 
and all J2 > 0.6| Ji|. Only close to the classical CAE/EM 
boundary at J2 = 0.5|Ji| does a two-magnon instability 
(with associated nematic order) prevail. These results 
are in complete agreement with previous exact analytic 
calculations for two-magnon bound states in the thermo- 
dynamic limit, and numerical exact diagonalizations of 
larger clusters'^. 

Using the two values (j) — 4'± ^^i the frustration an- 
gle together with the experimental energy scale Jc de- 
termined from zero-field susceptibility and heat capacity 
measurements2i2i2ii2£, we can extract the expected values 
for the two different saturation fields -ffgat = H± from the 
right-hand side plot of Fig. [3] Assuming a value g = 2 for 
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FIG. 3: Left: Scaling plot of the difference between the one- and two-magnon instability fields, as defined by Eq. (Illf) . for 
cluster sizes N — 16, 20 and 24 sites. The different symbols denote different positions in the phase diagram, corresponding 
to the values of listed on the right-hand side. Note that the size of the field difference is non-monotonous in (j>- Right; 
One-magnon (solid line) and two-magnon (dashed-line) instability fields for the 24-site cluster as a function of the frustration 
angle. In this and in subsequent plots of (^i-dependent quantities, the thin vertical lines denote the classical phase boundaries 
of the J1-J2 model^. 
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the average gyromagnetic ratio, we arrive at field values 
between 13 and 24 T, low enough to be reached experi- 
mentally. In Fig. [4] we have plotted the predicted values 
for Hsut as a function of the frustration angles (j)± for the 
known compounds (Fig.[T]), using the corresponding val- 
ues of Jc from zero-field measurements. Hsut can be de- 
termined for example by a magnetization measurement at 
sufficiently large fields. For the P04-based compounds, 
the saturation field together with the zero-field data for 
the susceptibility and the heat capacity would provide a 
direct way to determine the exchange constants Ji and 
J2 individually and hence the region of the phase dia- 
gram to which the compound belongs, without the need 
to measure the magnetic ordering vector directly. 



V 



FIG. 4: Predicted values for the saturation fields of the ex- 
perimentally known compounds. <!>+ labels the frustration 
angle corresponding to the collinear phase, and denotes 
the frustration angle for the Neel phase. The values for (j)± 
are determined from zero-field susceptibility and heat capac- 
ity measurementa^'^'^. 



B. Magnetization and susceptibility 



We have calculated the magnetization m(T, H) and the 
magnetic susceptibility x(T, iJ) = Np^^Q{dm{T,H)/dH) 
by evaluating the following thermodynamic traces: 



m{T,H) 



Jr 



1 Jr. 



N keT 



-{STY 



(13) 
(14) 



where we have explicitly included a factor 1 /N to account In the definition of the susceptibility, we also include the 
for the volume dependence of these extensive quantities. 
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magnetic permeability /io and the Avogadro number A^a- 
In order to make x(T, if) a dimensionless quantity, we 
need to multiply it with the characteristic energy scale 

At T = 0, the magnetization m{T = 0, H) of any 
finite size system evolves as a series of discrete steps. 
For a generic AF system with a singlet ground state, 
m{T = 0, H) takes on all possible (integer) spin values 
as a function of if, up to the field He at which the sys- 
tem saturates. Generally, in the thermodynamic limit, 
m(T = 0, ii < ffc) is a smooth curve, and singular fea- 
tures occur only where there is a magnetic phase transi- 
tion. However in the case of the spin-1/2 J1-J2 model, 
a "step" at exactly half the saturation magnetization 
m{T = 0, H) = 1/2 survives in the thermodynamic limit 
for J2 ~ </i/2, i.e. where a nonmagnetic ground state 
separates NAF and CAF order. This half-magnetization 
"plateau" is believed to be associated with the formation 



of localized magnon excitationssl. 

Temperature acts to smear jumps in magnetization. 
For the small clusters which we consider, the step- 
like behavior in m(T, H) has already disappeared for 
— 0.2 Jc. At the same time, all trace os the half- 
magnetization plateau is also lost. The magnetic sus- 
ceptibility therefore shows a smooth and nearly constant 
field dependence, see Fig. [S] It drops to zero upon reach- 
ing the saturation field. Only at the borders of the fer- 
romagnetic regime for h = do anomalies related to 
spontaneous magnetization appear. 



C. Entropy and heat capacity 

The entropy and heat capacity are defined through 



-^S(T,H) = l(l„Z(T,H) + J^(K(H))), (15) 

I 



again using the definition for the thermal averages given 
in Eqs. (|6l7p . Z and TC are the partition function and 
the Hamiltonian of Eq. ^ , respectively. The right-hand 
side of Fig.[5]shows a contour plot of the entropy S{T, H) 
at fixed temperature ksT = 0.2 Jc as a function of the 
frustration angle </> and the magnetic field H. In the 
ordered phases, the entropy is a smooth and almost con- 
stant function of the magnetic field, dropping to for 
fields higher than the saturation field iJsat- Character- 
istic anomalies can be observed at both edges of the 
coUinear phase, where the entropy crosses a broad max- 
imum as a function of field before again vanishing when 
crossing the saturation field. 

Figure [HI shows contour plots of the entropy (left) and 
heat capacity (right) as a function of magnetic field and 
temperature for a fixed value oi (f) — 0.747 tt. We have 
chosen this particular frustration angle for the plots be- 
cause it is believed to belong to = 0+ for the compound 
SrZnVO(P04)2- The wiggly contour lines at low temper- 
atures k-eT <C Jc are finite-size effects, where each tem- 
perature minimum corresponds to a Zeeman level cross- 
ing at T = as discussed in Section HV Al 

The lines of constant entropy (Fig. [51 left) are almost 
field-independent or even have a slightly negative slope as 
a function of field for low temperatures T <^ Jc/^b and 
fields H <SC iJsat- This implies that a sample cools down 
slightly when increasing H. The behavior of a paramag- 
net is opposite: Here, isentropic lines are straight lines 
crossing the origin, and a sample always heats up when 



increasing the applied field. Of course, for high enough 
fields and temperatures, the behavior of the entropy of 
the J1-J2 model is the same as that of a paramagnet. At 
the saturation field, which is iigat ~ 1-64 Jc/((7/xb) for 
(jy/n — 0.747, the temperature reaches a minimum when 
adiabatically changing the field at low temperatures, and 
rises steeply when increasing the field to higher values 
H > iisat- In this area of the phase diagram, a J1-J2 
compound is a good system for magnetic cooling, espe- 
cially in view of the low values for the saturation fields 
(in Tesla) for the experimentally known compounds (see 
Fig. [Hand its discussion above). 

The heat capacity Cv{T,H) is characterized by two 
maxima as shown in the right panel of Fig. [51 One max- 
imum occurs at if = and fceT « 0.5 Jc, which is the 
broad anomaly occurring at the crossover to the 1/T^- 
temperature dependence for high temperatures: We have 

This maximum has already been discussed in Ref. 2. A 
second maximum can be observed at high temperatures 
T ^ Jc/fce and magnetic fields H ^ iisat- Here, a field- 
induced gap opens, leading to a Schottky-type anomaly 
of an effective two-level system, see also Sec. IV Al and 
Fig. [m (center). 
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FIG. 5: (Color) Contour plots of the magnetic susceptibility x{T, H) (left) and the entropy S{T, H) (right) at a fixed temperature 
T = 0.2 Jc/fefi as a function of the frustration angle (j) and the magnetic field H. The plot was made using a 24-site cluster on 
a grid of 200 x 300 data points. 




2 3 
gy"BH/Jc 




2 3 



FIG. 6: Contour plots of the entropy S{T, H) (left) and the heat capacity Cv{T, H) (right) at a fixed frustration angle (p = 0.74 tt 
as a function of the magnetic field H and temperature T for a cluster of 24 sites. For the entropy plot, the contour line starting 
aX H = Q and k^T « 0.15 Jc corresponds to 5* = 0.05 A'^Afcs, the highest contour line starting at -ff = and k^T « 1.15 Jc 
corresponds to S = 0.6 A'aAib. In the plot of the heat capacity on the right-hand side, the lowest contour starting at iJ" = and 
fee ~ 0.1 Jc has a value of Cv ~ 0.05 NAks, while the highest contour starts at _ff = 0, ksT « 0.4 Jc and has Cv ~ 0.45 NAks- 



D. The magnetocaloric effect 

For the numerical calculation of the magnetocaloric 
effect we express Eq. ([T|) as the cumulant 



dT 
M 



{TP) - {nf 



(18) 



and normalize the results to the magnetocaloric effect 
of a paramagnet. The left-hand side of Fig. [7] shows a 
contour plot of Fmc — ^inc/ {T/H) as a function of the 
frustration angle (j> and the magnetic field H . For small 
fields H <^ Hsat, ^mc/ (T/H) is small, nearly zero or even 
slightly negative, apart from finite-size effects showing up 
in particular in the coUinear phase. It is only at the sat- 



uration field where T^^/ {T/H) develops a large anomaly 
peaked slightly above i?sat (compare with the right plot 
of Fig. [3]). For magnetic fields H ^ i?sat, we eventually 
reach T^,,/{T/H) ^ 1. 

Apart from a factor T, the magnetocaloric effect is 
given by the ratio of two quantities, see Eq. ([T]): (a) 
In the numerator, we have {dM/dT)H, or, equivalently, 
{dS/dH)T- The entropy S{T,H) at constant tempera- 
ture is plotted on the right-hand side of Fig. [5l Its field 
dependence, corresponding to the density of contour lines 
in the plot, is weak apart from the nonmagnetic regions 
at the edges of the collinear phase, (b) The denomina- 
tor is the heat capacity CviT,H). Here, small values 
give rise to a large magnetocaloric effect. Fig. [7] (left) 
holds a plot of the heat capacity at constant temperature 
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FIG. 7: (Color) Contour plots of the heat capacity (left) and the normalized magnetocaloric effect T^c/ {T/ H) (right) for the 
24-site cluster at a fixed temperature T = 0.2 Jc/fce as a function of the frustration angle and the magnetic field H. 




FIG. 8: Specific heat at ksT/Jc = 0.2 and constant, tj)- 
independent field gfj,^H/Jc ~ 2.5 as function of frustration 
angle. The double peak structure at the classical NAF/CAF 
boundary corresponds two the two ridges in the contour plot 

of Fig. m 



T = 0.2 Jc/fce as a function of the frustration angle (j) and 
the magnetic field H. The heat capacity is large in the 
disordered regions, reflecting the high number of quasi- 
degenerate states. Around J2/J1 = 1/2 {(P/tt « 0.148), a 
two-peak structure evolves when increasing the field, see 
also Fig [51 Due to the smallness of the saturation field, 
we currently cannot say whether such a structure also ex- 
ists at the "mirrored" (J2 —J2) position in the phase 
diagram at J2 / Ji =—1/2. When reaching the saturation 
field, the heat capacity drops and eventually vanishes. 

Taken together, it appears naturally that the magne- 
tocaloric effect is peaked around the saturation field. The 
drop in magnitude inside the nonmagnetic regions can be 
understood, too, as a consequence of their large specific 
heat. And since the entropy rises inside these regions 



when turning on the magnetic field, the magnetocaloric 
effect must be negative, indicating a cooling of a sample 
before reaching the entropy maximum. We note that we 
also observe a change of sign in {dS/dH)T as a function 
of field in the magnetically ordered regions, and return 
to this point in the context of spin wave theory below. 

On the right-hand side of Fig. [51 we have plotted the 
values of r,iic(r, H) at the saturation field as a function 
of the frustration angle again for T = 0.2 Jc/fce- The 
open circles denote the absolute values, while the filled 
circles denote the values relative to a paramagnet. In ac- 
cordance with the discussion in the previous paragraph, 
the absolute values of T^c{T,H) in the ordered phases 
are larger than in the nonmagnetic regions. The devi- 
ation from the average value is less than a factor two. 
In contrast, the normalization to the magnetocaloric ef- 
fect of the paramagnet (filled circles in the right panel 
of Fig. ^ introduces a strong influence of the saturation 
field, compare the right plot in Fig. [31 Therefore, the 
highest enhancement of Fmc with respect to a paramag- 
net occurs deep inside the magnetically ordered phases, 
where the saturation field reaches its maximum values. 

The heat capacity Cy(T, ffgat) as a function of the sat- 
uration field at constant temperature T = 0.2 Jc/^b is 
plotted on the left side of Fig. [51 It is strongly enhanced 
in the nonmagnetic regions, while roughly constant as a 
function of the frustration angle (p in the magnetically 
ordered phases, giving rise to the comparatively weak (p 
dependence of Fmc- As discussed above, the enhance- 
ment inside the disordered regions is responsible for the 
suppression of the magnetocaloric effect. 

To clarify the field dependence of V^c{TtH) further. 
Fig. [TUl holds a comparison of the three relevant quan- 
tities T^,{T,H)/{T/H) = {dT/dH)/{T/H), Cv{T,H), 
and S{T, H) as a function of the magnetic field H for 
a fixed frustration angle (p = 0.747 tt and a temperature 
T = 0.2 Jc/fce: Disregarding possible finite-size effects, 
the entropy (solid line, left scale) and the heat capacity 
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FIG. 9: In this figure tlie field is kept at tfie saturation field for every (j>, i- e. Cy and Fmc are plotted along the curve in Fig. |3] 
(right). Left: Value of the heat capacity Cv{T,H) at the saturation field as a function of the frustration angle </!> for a fixed 
temperature T = 0.2 Jc/fee- Right: Values of the magnetocaloric effect Fj^ciT, H) = (dT/dH)s at the saturation field as a 
function of the frustration angle (j> using a fixed temperature T — 0.2 Jc/fca. The open circles denote the absolute value of Fmc, 
their scale is given at the left ordinate. The filled circles denote the normalized values Vyac/{T/ H), indicating the enhancement 
relative to a paramagnet. Scale is on the right ordinate. Both plots were made using a 24-site cluster. 



(dotted line, left scale) are slowly varying functions of H , 
see also Sec. IV A[ dropping sharply above the saturation 
field Hsa,t ~ 1-64 Jc/((7/xb). Taken together, this leads 
to a pronounced maximum of {dT / dH) / {T / H) slightly 
above the saturation field. Otherwise, [dT / dH)/ {T / H) 
is small and negative for fields H ^ Hgat (because 
{dS/dH)T > in this field range) and approaches 1 for 



V. APPROXIMATE ANALYTICAL 
TREATMENTS OF THE MODEL 

To better understand the exact numerical results for 
finite clusters it is useful to have approximate analyt- 
ical results available for comparison. We consider two 
approaches: Firstly a mean field treatment which pro- 
vides a reference point for the global behavior of entropy 
and specific heat in the ordered phase, and the magne- 
tocaloric effect above the saturation field. Secondly we 
use a linear spin-wave (LSW) approximation to investi- 
gate the anomalous enhancement of the magnetocaloric 
effect around the saturation field, which turns out to 
be due to a softening of spin excitations at character- 
istic wave vectors. In this approximation, in contrast to 
mean field theory, the MCE below the saturation field 
is nonzero. The spin wave approximation also allows to 
study subtle effects for subcritical fields which lead to a 
sign change of the MCE. 

The existence of long range order at finite tempera- 
tures, i. e. a non-vanishing transition temperature Tc as- 
sociated with magnetic order, is implicit in both treat- 
ments. In reality, for the layered vanadates which we wish 
to describe, Tc will be determined by low energy scales 



which are not present in our model, notably the interlayer 
magnetic exchange Jj^ and magnetic anisotropy 5J. 



Formally, we cannot break a continuous symmetry such 
as spin rotation at any finite temperature in 2D, and to 
be "correct" we should generalize our model to higher 
dimension and finite anisotropy. However these details 
make little qualitative (or quantitative) difference for a 
wide range of temperatures SJ,J± <^ T <C so we 
suppress them below. Furthermore, in these calculations 
we neglect the consequence of interactions between spin 
waves (see e.gi^) and the break-down of magnetic order 
at zero temperature on the borders of the CAF phase^iS^. 
These effects can be expected to modify the details of 
critical behavior as a function of magnetic field, but not 
its broad features, and are left for future investigation. 



As discussed in the previous section, the size of the 
magnetocaloric cooling rate in Eq. ([1]) is determined by 
the ratio of the rate of change of entropy with field 
{dS I dH)rp y and its rate of change with temperature 
{dS/dT)^ y = Cv/T. On approaching the critical field 
both quantities tend to increase sharply, and the result- 
ing increase in T^c{H) is a tradeoff between them. It 
is not immediately obvious for which frustration angle 
(j) the enhancement in Fnic(^^c,<^) should be largest. At 
modest temperatures, the simple spin wave approxima- 
tion described below gives considerable insight into this 
question, and the critical anomalies of the MCE around 
the saturation field. 
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qHbH/Jc H/Hsat 

FIG. 10: Left: Entropy S{T,H) (solid line, left scale), heat capacity Cv(T,H) (dashed line, left scale), and magnetocaloric 
effect (dT/dH) / {T/ H) (solid line, right scale) as functions of the magnetic field H at constant temperature T = 0.2 Jc/fe for 
a frustration angle (j) = 0.74 tt. Data were generated using a 24-site cluster. Right: Normalized MCE for various temperatures. 
The anomaly at He is suppressed with increasing T, see also Fig. [TS] At the lowest temperature a negative MCE is possible. 



A. Calculation of mean field order parameters and 
thermodynamics 

In this section the magnetothermal properties will 
be investigated in mean field approximation to have 
a reference for the spin wave and numerical exact- 
diagonalization methods. The results of the former are, 
however, not expected to give a realistic description of the 
MCE. For a unified treatment of AF phases it is advisable 
to use a four-sublattice description {a,f3 — A,B,C,D) 
with each sublattice having N/4 sites for both NAF and 
CAF. Since we consider only isotropic exchange we may 
assume without loss of generality that the field is per- 
pendicular to the (xy) plane of the square lattice, i.e., 
h = hz. 

In this and the following subsection we refer all exten- 
sive quantities like entropy, specific heat etc. to a single 
site for convenience. The exchange field h"'^ and total 
molecular field h due to Eq. Q is then given by 

h„ = h + h"'^ (19) 

where the exchange constants J^^ are defined per bond. 
The components of the exchange field h'^^ and /i]^ which 
are parallel and perpendicular to the field direction z are 
related to the respective spin expectation values (S'[|) and 
{S±_) via the equations 

Iftf = -a||(5||); i/i- = a^(5x> (20) 
where the prefactors for the AF and the FM or fully 



polarized phases {h > he for any 0) are given by 

NAF: ail = |(Ji + J2); - |(Ji - J2) 

CAF: an = |(Ji + J2); «± = |>^2 (21) 

FM: an = |(Ji + Jj); ax - 

Then the mean field approximation of the Hamiltonian 
in Eq. ^ may be written as 

il,„f = ^[-(h + h-)S:, + ih-(S„)] (22) 

This also defines selfconsistently the mean field aver- 
ages via (A) - Tr{Aexp(-/3i/,M)}/Tr{exp(-/3H,M)} 
with f3 = l/lk-oT). From the above equation the corre- 
sponding total mean field internal energy per site J/mf — 
{\/N){Hraf) is obtained as 

[/„,f(T, H) = --J2{\,+ ih-)(S„) (23) 

Explicitly, using Eq. pT|) one obtains 

NAF, CAF: C/^f = -h{S\^ + a\\{S\\f - a^{SA_f 

FM: C/^f = -h{S{) + a\\{S\\f (24) 

To calculate thermodynamic quantities the expecta- 
tion values {S\\) and {S\_) and their temperature deriva- 
tives have to Idc obtained selfconsistently. This is done 
by diagonalizing i/mf which leads to local eigenstates 
|±) = "±1 T ) + ^^±1 i ) where | t ), I i ) are the degen- 
erate free S = 1/2 states. The |±) states have energies 
E± = Ec + e± given by 

= ax(5x)2-a||(5||)2 

with h\\^h + hf; hi_ ^ h°l. (25) 
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FIG. 11: Field dependence of staggered order parameter {Sx) and uniform magnetization {Sz) for two different temperatures 

(left and center panel). The total moment {S) — {{Sx)^ + {Sz)'^)^ is field independent in the ordered regime. The field 
dependence of the canting angle Oc/2 (counted from the field direction) is also shown. Right panel: Susceptibility x for three 
different temperatures. Its value in the ordered regime is T- independent. In all cases the frustration angle is tf) — O.TAn 
corresponding to the CAF choice. 



Defining 2eo = A — e+ — e_ as the splitting due to the 
molecular field h their coefficients are then obtained as 



u± 



[eo(2eo±/i||)]5 



v± 



Using these coefficients and the difference of thermal 
occupation numbers p- — p+ = tanh j/JA of eigenstates 
|±) one finally obtains selfconsistent equations for the 
spin expectation values: 



{S±) = i^tanhi/3A, 



Ih 



1 



= ^tanh-^A 
^ 2 A 2 



(27) 



The selfconsistency is implied via the molecular field ex- 
pressions 



= h — 2a[| {S\\ 
h± = 2aj_{S±) 



(28) 



The canting angle Oc/2 of magnetic moments is obtained 
from minimizing U^nf m Eq. (j24[) The angle is counted 
from the field- or c-direction and given by 



or 



tan(-^) 



.^c, h 1 



h 

hr 



(29) 
(30) 



where hc=2(S)(a|| -l-aj.). Inserting this into Eqs. (|27I28|) 
leads to the simple and general result 

(5) = itanhi/3A with /^^2{S)\ax_\ (31) 

This means that in the ordered phase the molecular field 
splitting A of spins is field independent up to and 
hence the total moment {S) is also field independent, i. e. 
the moment can only be rotated by the field as long as 
the transverse staggered order exists. This fact has strik- 
ing consequences for the thermodynamic quantities below 
/ic- For the thermodynamics we also need the tempera- 
ture derivatives 9(S',)/ar = -fcB/32(i9(5',)/a/3) {i =||,-L). 
They are obtained from Eq. (|27p in a straightforward but 
lengthy calculation and the resulting explicit expressions 
are given in appendix A. 

The mean field solution for a CAF value of (p = 0.74 tt 
is shown in Fig. [TT] On the left panel the decrease of 
the staggered OP {S^) with increasing field and the con- 
comitant increase of the uniform moment {Sz) are shown 
for small temperature. The total moment {S) is constant 
as predicted and practically equal to S=l/2 in the whole 
field range. This confirms that the moment is simply ro- 
tated (canted) by the field without changing its size. The 
relevant canting angle 9c/2 is also shown in the figure. 
For moderate temperatures (center panel) the zero field 
value of {Sx) is already reduced somewhat. As required 
by Eq. (PT|) the field still only reorients the moment, i. e., 
{S) is a constant less than 1/2 for fields h < he- Fi- 
nally for h > he when the moment is aligned with the 
field {9^/2 = 0), the total moment (S) = (S^) will be 
polarized, i. e. increases with field until it reaches the 
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FIG. 12: Thermodynamic properties in mean field approximation for = 0.747r (CAF) as function of applied field for three 
subcritical temperatures fcaT/Jc — 0.3 (a), 0.5 (b) and 0.7 (c). Left: Entropy dependence on h. Below the saturation field 
he it is field independent because (5*) (see Fig. Illjl and also A is field independent in the ordered regime. Center: specific 
heat dependence on h. For larger temperatures when he is sufficiently suppressed a Schottky peak evolves above he. Right: 
Magnetocaloric enhancement factor. In MF approximation almost no enhancement is visible above hc{T). Below the critical 
field the entropy is field independent and hence Fmc drops to zero suddenly according to Eq. ([TJ. 



asymptotic value of 5 = 1/2. The right panel shows the 
susceptibility Xmf for various temperatures. It is constant 
and T-independent in the ordered regime. These results 
are qualitatively unchanged for different angles 0. 

The desired thermodynamic quantities may now be 
conveniently obtained from the free energy F,„f and in- 
ternal energy Umf fEg. lM)) of the S* = 1/2 system split by 
the molecular field by an energy A(/i||, hj_). The former 
is given by 

F,„f (T, H)=E,--^ ln[2 cosh i (/3A)] (32) 



per site and the entropy 5mf — —{dFnii/dT), specific 
heat Cy^ = {dUmf/dT) and susceptibility per site are 
then obtained as 



Xmf 



fee [ln(2 cosh i/3A) 

fcB/3'p-2a||(5|i»(^ll 

{9^Ji^f^{S\\)/^h 



i/3Atanhi/?A] 



where (<S'||)', {S\_)' are given in Eq. (jAip of Appendix A. 
For uncoupled spins ay — = Q and {S\\)' — 

■j/icosh"^ which leads to the Schottky specific heat 
of the two level system. For the magnetocaloric cooling 
rate Fmc we need in addition the temperature gradient of 
the magnetization Wmf = g^J'B{S\\) which is simply given 
by 



dT 



-te(5MB)/3'(5||)' 



(33) 



Then the mean field expression for the cooling rate Fmc 
may be obtained from the definition in Eq. ([1]) using the 
expressions for Cy in Eq. (|33p and the temperature gra- 
dient in Eq. ([33]) . Thus the solution of the selfconsis- 
tent Eqs. (P7|) for (5*11) and and their temperature 

gradients and (5^)' in Eq. (|A1|) provide all the 

necessary input for obtaining the magnetocaloric quanti- 
ties S'mf(r, H), CJ^f (T, H) and T^^{T, H) from the above 
equations. 



In Fig. [13 we show the field dependence of mean field 
entropy, specific heat and cooling rate as function of field 
for various temperatures. Obviously S and Cy (left and 
center panel) are field independent below caused by 
the fact that A is constant according to Eq. ([5T|) . Con- 
sequently, since the cooling rate is proportional to the 
field-gradient of S (Eq. [1]) , it suddenly drops to zero be- 
low the saturation field he as seen in the right panel of 
Fig. [T21 In addition it shows that above he the cool- 
ing rate is only slightly enhanced from the paramagnetic 
value. Although the mean field results are far from real- 
istic, the two main aspects, field independence of S and 
Cv below he and steplike anomaly in Fmc at h^ still leave 
their signature in the more advanced spin wave and nu- 
merical treatment to be discussed below. 
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B. The magnetocaloric effect in the linear spin 
wave approximation 

1. Statistical mechanics and general expressions 

In mean field approximation the elementary excita- 
tions are local dispersionless spin flips whose energy is 
determined by the molecular field. This is far from real- 
ity especially close to the saturation field when the order 
parameter breaks down, associated with a softening of 
the spin excitations at some k-point or even line in the 
Brillouin zone (BZ). The linear spin wave approxima- 
tion takes this effect into account and gives a much more 
realistic description of the MCE. Since spin wave interac- 
tions are left out in our approach, however, the singular 
behavior around may be overestimated. 

To calculate the magnetocaloric effect of Eq. ([T]) in spin 
wave approximation we start from the partition function 



Tr[e" 



'■H/ksT-i 



(34) 



where the Hamiltonian Ti. is expanded in spin wave coor- 
dinates using the Holstcin-Primakoff approximation 



Ak 

+ OiEo/S^) (35) 



Here Eq = Umf{T — 0) is the classical (mean field) 
ground state energy per spin, ex{h,\i) the spin wave dis- 
persion in applied field h, and the sum over k runs over 
the appropriate magnetic BZ, while A counts the differ- 
ent spin wave branches within that BZ. The operator a^^. 
creates magnons with commutation relations 



'Ak' "A'k 



In addition to the classical (mean field) ground state en- 
ergy per spin Ea there is zero point energy contribution 



^^p = ^E[^A(/i,k)-^(/i,k)] 



(36) 



Ak 



where A{h, k) is the on-sublattice coupling between spins, 
defined below. For CAE and NAF phases the E^p is 
always negative; in the EM, where the ground state 
and spin waves are eigenstates with a single dispersion 
e{h, k) = A{h, k), Ezp vanishes identically. 

The partition function is essentially that of set of in- 
dependent simple harmonic oscillators 



Z = g-[Eo{h) + E^^ih)]/ik^T) 



xn[i 

Ak 



-eA(/i,k)/(feBT) 



(37) 



From this we find the internal and free energy per site. 



using nB(e,T) = [e'^/C^^T) - i] i for the Bose factor: 



N ■ 

A,k 

F = -— fceTlnZ 
N 

= Eo{h)+E,p{h) 
+ 4 fceT^lnJ 



(38) 



N 



1 



Ak 



The entropy per site S 



-{dF/dT) follows directly 



A,k 



_ctnh^^-lnsinh^4^ 



(39) 



We can also find the uniform magnetization 

_ dE^jh) dE,p{h) 
^ dh dh 



Ak 



dh 



(40) 



In the EM phase, where e{h, k) = uj^ + h and there is no 
zero-point term, this simply reduces to 



m — Too 



^^nB(wfc + /i,T) 

Ak 



(41) 



with Too = dEo{h)/dh. 

Quite generally, we can calculate the MCE as the ratio 
in Eq. ([1]) where the magnetization gradient is given by 



dm 
'dT 



1 

'n 



4(fcBT)2sinh^ 



2fcBT 



(42) 



and the specific heat Cy = T{dS/dT) ^ [dU/dT) by 



Cv _ kB MT-^ 



ex{h,k) 


2 


4(fcBr)3sinh2 


2kuT 



(43) 



From Eqs. (gSlllSl), Fmc is obtained using Eq. This 
reduces the problem to one of evaluating two-dimensional 
integrals on k for the appropriate spin wave dispersion 

We note that in the special non-interacting case 
e(/i, k) = /i, these expressions reduce to those for an 
ideal quantum paramagnet, with the associated magne- 
tocaloric effect: 



dT 
dH 



T 
H 



(44) 



Incidentally this is a general property of any system for 
which the partition function Z depends only on H/T, i. e. 
F = ~kBT\nZ{H/T). 
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FIG. 13: (Color) Spin wave dispersions in the NAF and CAF phase. Only the e+(/i,k) modes are shown. The corresponding 
e_(/i,k) modes are obtained by translation with (tt, tt) or (vr, 0) in the NAF and CAF case respectively. Left panel: canted 
NAF dispersion e+(/i,k) for J2 — and a canting angle of 9c = 7n/8. Note that e+(/i,k) is gapless at (tt, tt) and gapped at 
(0,0) while the opposite holds for e_(/!,,k). Right panel: canted CAF dispersion e+(/i,k) for Ji = 1, J2 = 1, S = 1/2 and a 
canting angle of 6c = 7n/8. In this case e+{h,'k) is gapless at (7r,0) and gapped at (0,0) and vice versa for e_(/i,k). 



2. Spin wave dispersion in the ferromagnet/ saturated 
paramagnet 

Since spin wave theory assumes broken spin rota- 
tion symmetry, we can treat the spontaneously polarized 
FM phase for h — and the saturated paramagnet for 
h > hc{(j>) on an equal footing. Expanding about the 
maximally polarized state to 0{1/S) we find 

H = NEo + "^eFM{h,k)alaj^ + 0{l/S^) (45) 



where 



Eo = 2( Ji + J2)S^ - hS 



(46) 



is the classical ground state energy per spin, which is 
equal to Urai{T — 0) in Eq. ([M)) . The spin wave disper- 
sion has a single branch given by: 

e^mih, k) - - 4Ji5'[l - 7(k)] - A.hS[\ - 7(k)] (47) 

where 

7(k) = ^(cosfca; -f cosfcj^) and (48) 



7(k) = cos kx cos ky 

Wc then have simply 

deYM{h,'k) _ 
dh 



(49) 



(50) 



i. e. a rigid shift of the entire dispersion with change in 
magnetic field. 

The dispersion will have a single (parabolic) mini- 
mum at k = (0,0) in the FM phase (i.e. for Ji < 0, 
J2 < |>^i|/2). However in the saturated paramagnetic 
state, above the critical field hc{4>,Jc), the minimum of 
the dispersion will be at k = (tt, tt) where there is a NAF 
ground state, and at k = (tt, 0) (and symmetry points) 
where there is a CAF ground state. At the classical 
critical point Ji — 2J2 > separating NAF from the 
CAF, there are line zeros around the zone boundary 
kx = ±7r and ky = ±7r. At the classical critical point 
— Ji — 2J2 > separating FM from the CAF, there are 
line zeros for kx — and ky = 0. Note that these line 
zeros connect the different minima of the dispersions 
between which this state must interpolate. 



3. Spin wave dispersion in the canted NAF 



Expanding about a canted NAF with ordering vector 
(7r,7r) and canting angle 9c we find 



n^NEo+J2 [^(^^ k) (a^a^ + 6[& J + B{h, k) (aj, foL^ + ^k^-k) + C{h, k) (aj^fo^ + &j,a J ] + OiEo/S^) (51) 

k 
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The classical ground state energy per site is given by 

Eo = 2JiS^cos9c + - hS cos{0j2) (52) 

which is identical to U^iiT — 0) of Eq. (IM)) . Minimizing 
this energy fixes the canting angle (Eq. (1301) 1 

(53) 



where {dc/2) is measured relative to the magnetic field 
direction z (i. e. the FM has 8c = 0, the NAF 9^ = n/2). 
Exactly the same expression follows from the requirement 
that the spin wave expansion contains no terms linear in 
bosons. Eliminating the magnetic field through Eq. ([53]) . 
we find: 

A{h,k)±C{h,-k) = 4JiS' [l±cos2(6ic/2)7(k)] 

- AJ2S [1 - 7(k)] (54) 
B(/i,k) = -4JiS'7(k)sin2(6'c/2) 

These expressions still depend on the applied magnetic 
field through the canting angle Of.{h). The bilinear form 
of Eq. (|5ip (Appendix B) may be diagonalized by a Bo- 
goliubov transformation to give 

n = NEn+NE,p+ J2 «A(/i,k)alka^k 

k,A=± 

+ 0(1/8^) (55) 

where 



e±(/i,k) = ViMh, k) ± C{h, k)]2 - B{h, k)2 (56) 

The two-fold degeneracy of the spin wave dispersion 
of the NAF is lifted by the applied magnetic field. 
In the (physical) magnetic Brillouin zone centered on 
k — (tt, tt), the Goldstone mode is e+(/i, k) = 0, while 
e_ (h, \i) = h has a finite gap. Since spin waves are not 
eigenstates, there is now a zero-point energy term in the 
energy Ezp. 

In order to calculate the rate of change of magnetiza- 
tion with magnetic field in Eq. (I42p . we also need the 
field derivative of the dispersion. Considering explicitly 
e+(/i,k), we obtain 



de+{h,k) _ 1 2Ck 
dh ~ 7+~h' 



(57) 



As far as the periodicity of dispersions is concerned, note 
that translation by Q=(7r,7r) leads to 

A{h,k + Q) + C{h,k + Q) = A{h,k} -C{h,k) 

B(/i,k + Q) = -B{h,k) (58) 

therefore the two e± modes are simply interchanged by 
translation through the NAF ordering vector Q=(7r, tt). 
This means that in simple thermodynamic averages one 
can work with a single mode in the full square lattice 
BZ — e. g. e+{h, k) — rather than with the two (physically 
distinct) modes in the smaller magnetic BZ. The e^{h, k) 
spin wave dispersion for NAF in the full paramagnetic BZ 
is shown in Fig. [13] (left panel). 



4. Spin wave dispersion in the canted CAF 

The ground state energy per spin is now given by 

Eo = JiS^[l + cos9c]+2J2S^cosec 

-hScos{6c/2) (59) 

which again is equal to Unii{T — 0) in Eq. ((32)) . Mini- 
mizing this energy leads to a canting angle (cf. Eq. [5D|) 



— = cos 
2 



h 



4:JlS + 8J2S 



(60) 



Once again we obtain spin waves with a dispersion of 
the form Eq. ([5^)1 . After elimination of the field using 
Eq. ([50]) . the coefficients of the spin wave expansion are 
given by 

A{h,k) = 25[2J2 + Jicosfcj,] (61) 
B{h,k) = -2S[Ji + 2J2 cos ky] cos k^sin^{0c/2) 
C{h,k) = 25[Ji + 2J2Cosfcj^]cosfc:^:Cos^(6'c/2) 



These coefficients once again satisfy the relation Eq. (|58|) 
with Q=(7r, 0), i. e. the e± modes are interchanged under 
translation through the magnetic ordering vector. There- 
fore simple averages can again be calculated for a single 
mode in the full paramagnetic (square lattice) BZ. As in 
the NAF case, the field gradient of spin wave energies is 
given by Eq. (|57)l. The e+(/i, k) spin wave dispersion for 
the CAF in the full paramagnetic BZ is shown in Fig. [T51 
(right panel). 



VI. DISCUSSION OF THE ANALYTICAL 
RESULTS AND COMPARISON WITH THE 
NUMERICAL FINDINGS 

The typical field dependence of entropy, specific heat 
and MCE in the linear spin wave approximation as cal- 
culated from Eqs. ([H EH [H 031) are shown in Fig. [H 
Entropy, specific heat and MCE are all smooth functions 
of magnetic field except at the critical field at which 
there is a (2nd-order) phase transition between the para- 
magnet and canted Neel phases. 

The most striking feature of these predictions is the 
double spike in the MCE at h^. This is a generic feature 
of a 2nd order phase transition between paramagnetic 
and ordered phases in applied field^'^ and has previously 
been seen in Monte Carlo simulations of the classical 
Heisenberg model^S,. This sudden and sharp change in 
sign of the MCE can easily be understood in terms of 
the contours of fixed entropy (adiabats), discussed be- 
low. It is accompanied by closely related cusps in the 
entropy and heat capacity, peaked at he. 

In general, as we would expect, entropy and specific 
heat have a much weaker field dependence in the ordered 
phase below he than in the disordered phase above it. 
The entropy of the (gapped) paramagnetic phase falls 
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FIG. 14: Left panel: Entropy, specific heat and MCE for (f> in the CAE regime. The MCE has a maximum at ftmax and then 
changes sign at ho. At the saturation field /ic another sign change associated with a spike- like singularity appears due to the 
field gradient of the entropy. Right panel: Normalized MCE cooling rate Tmc/^mc ~ {H /T)Tjac for five different frustration 
angles in the NAE (0 = — 0.257r, 0) and CAP {<j) — 0.257r, O.Svr, O.TStt) regime and shown for moderate fields. The field is 
normalized to the saturation field hc{<j>) given in Fig. [3] (right). Low field maximum of the MCE at /imax(<^) and sign change at 
ho{(j}) are clearly seen to occur for all frustration angles. Pairwise equalities of the MCE are observed due to the symmetries of 
the spinwave spectrum with respect to (j>. 
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FIG. 15: Contour plots of the entropy S(T,H) in mean field approximation (left panel) and spin wave approximation (right 
panel) in steps of AS = 0.02A'AfcB and 0.05 A'a fee respectively. The former shows a temperature dependent critical field He 
but no structure of the entropy. The latter has a constant critical field but exhibits the cusp structure around He which is 
responsible for the sign change and peaks of the MCE in Fig. 1141 These contour plots should be compared with the results 
from FTLM in Fig. [6] (left panel) which qualitatively exhibit both features. 



rapidly in applied magnetic field, while the Neel phase 
responds to magnetic field by canting, at nearly constant 
entropy. As a result the typical (absolute) value of the 
MCE is much larger above he than below. These features 
of the LSW predictions are reminiscent of the mean field 
theory, as illustrated in Fig. [T^l 

For sufficiently large h ^ T, Jc we must (and do) re- 
cover the response of an isolated paramagnetic spin in 
either approximation. Where the LSW predictions differ 
from those of mean field theory is in the singular fea- 
tures at he, and in the presence of a finite MCE in the 
ordered phase. This is positive for h ^ 0, exhibits a shal- 
low maximum at a characteristic field h = /imax, changes 
sign at another characteristic field h = ho, before exhibit- 
ing the dramatic double spike at h = he — irrespective of 



which ordered phase is considered. The absolute values 
of these critical fields, and the form of the anomalies at 
he, do however depend on the structure of the low energy 
spin spectrum, and therefore on frustration through (/). 

The difference between mean field and spin wave re- 
sults becomes most obvious in a comparison contour plot 
of the entropy S{T, H) shown in Fig. 1151 The former 
(left panel) has a temperature dependent critical field 
which vanishes at the (mean field) transition tempera- 
ture. The entropy does not have any structure below he 
due to the field independent molecular field splitting A 
in Eq.( I3ip . In spin wave approximation (right panel) 
the critical field is temperature independent, but the en- 
tropy contours show a typical cusp structure around he 
with a maximum at ho further down which is caused by 
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FIG. 16: Thermodynamic quantities slightly above the critical field ^cCt/") as function of <j) for temperature k^T/Jc = 0.2. Here 
we used Ah = h — he = 10~^hc. Left: Specific heat shows enhancement at the NAF/CAF boundary. Right: Bare (dashed 
fine) and normalized (solid line) magnetocaloric effect as function of (p. Giant peak at the NAF/CAF boundary occurs due 
to gapless LSW modes along lines in the BZ. The overall behavior of the normalized Fmc/F^Jje — {Hc{(t)) /T)Tn-ic follows the 
f/i-dependence of the critical field (Fig. [Sj right)). 



the excitation of low energy spin waves (cf. Ref. [25ll26l ). 
According to the definition of Fmc in Eq. ([1]) this imme- 
diately translates into the sign change of r,nc at /iq and 
its enhancement around h^. 

The behavior of the (normahzed) MCE in the low to 
moderate field regime is presented on an enlarged scale 
in Fig. [M] (right panel) for typical values of (p in the 
NAF/CAF region. We notice a considerable variation 
of the characteristic fields ft.max(0) and /).o(0) with the 
frustration angle. The maximum MCE, rnic(ft.max) for 
ksT/Jc = 0.2 is of the order of ten percent of the para- 
magnetic value F^j.(/i,jiax)- Furthermore a symmetry in 
the 4> dependence is obvious: Firstly the MCE is invariant 
under reflection at the axis cj) = O.Stt or Ji = when both 
values of 4> lie in the CAF sector. This is obvious from 
the spin wave dispersion Eq. (|56l6ip in the CAF regime 
which is invariant under the simultaneous transforma- 
tion (Ji, J2) (- Ji, J2) and {K,ky) (fc^; -I-tt, ky + ir). 
Since the MCE is obtained by integration over the whole 
BZ, r„ic is unchanged under sign reversal Ji — > —Ji. Sec- 
ondly the MCE for (/) = and = O.Stt are equal, i. e., it 
is invariant under the replacement ( Ji, 0) — > (0, J2). 

Finite size effects prevent the characteristic fields 
hraayi{<i>) and /io(0) from being identified in FTLM cal- 
culations, as shown in Fig. [Tni Nevertheless a negative 
MCE at moderate fields is clearly compatible with the nu- 
merical results. In practice, for a finite size cluster, each 
of the ground state level crossings shown in Fig. [5] shows 
up as a separate "phase transition" in the FTLM re- 
sults for the MCE, with associated positive and negative 
spikes in Fig. [TOl For this reason, the sign of the MCE 
remains ambiguous. However for the saturated paramag- 
netic phase, where there is no further level crossing in the 



ground state, the sign of the MCE is correctly resolved 
and a pronounced enhancement is seen approaching hc{4>) 
from above. 

This is in good qualitative agreement with the pre- 
dictions of LSW theory, where it is clear from Fig. [T3] 
that the largest positive MCE is to be expected just 
above hc{4>). This maximum arises from the closing of the 
spin wave gap in the fully polarized phase when h h'^ . 
It occurs for fields slightly above /ic(</>) because temper- 
ature acts to "round" the sharp cusp in the entropy con- 
tours at the critical field — cf. Fig. [151 Needless to say, 
the low energy excitations responsible for the sharply di- 
verging peaks seen in LSW results are not accurately 
described by a cluster of 24 sites, and so it makes lit- 
tle sense to compare FTLM and LSW predictions at a 
quantitative level. 

While the structure of the MCE is generic to a (sec- 
ond order) phase transition, the details depend strongly 
on the amount of frustration present. If the frustra- 
tion angle is deeply within one of the ordered sectors, 
the softening of the spin waves at /i — > /i+ occurs at 
the wave vector of the low field AF structure (7r,7r) or 
(0,7r). However if the frustration angle approaches the 
transition regions CAF/NAF and CAF/FM, the soften- 
ing will occur along the whole line in the BZ connecting 
the wave vectors of the competing structures. This is 
reminiscent of, but less dramatic than, the situation in 
certain geometrically frustrated magnets where the spin 
gap closes simultaneously at he for an entire branch of 
excitations across the BZ, leading to the condensation 
of a macroscopic number of localized magnon modes^^. 
In the present case one should therefore expect a strong 
enhancement of Tync{h ^ /i^) for (p close to one of the 
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above boundary regions. 

The same holds for the specific heat. In Fig. [11] (left 
panel) we show the peak value Cv{h'^, <P) as function of 
frustration angle. Indeed the specific heat shows a strong 
enhancement for <j) ~ O.IStt (CAF/NAF) and ~ O.SStt 
(CAF/FM) of considerable width in (p. This is in good 
qualitative agreement with the FTLM results for finite 
clusters presented in Fig. [5] (left panel). The anomaly 
at the FM/NAF boundary on the other hand is much 
smaller. 

On the right panel of Fig. [11] the corresponding plot 
for the MCE is shown. The dashed line shows the bare 
MCE coefficient rmc(^ h'^) as function of (f>. It is 
almost constant except at the phase boundaries where 
again a large, but much narrower peak appears. This 
is not immediately obvious since C'v enters in the de- 
nominator of the expression for in Eq. ([1]). In fact 
on approaching ^/tt ~ 0.15 the MCE slightly decreases, 
only very close to the value when the spin wave dispersion 
softens along the line (tt, 0)-(7r, tt) and equivalent ones in 
the BZ a very sharp spike appears. This is due to the fact 
that the spin wave softening leads to a stronger increase 
of the magnetization gradient (Eq. ((12|)) sls compared to 
the increase of Cy from Eq. (|43)) . The full line in Fig. [TBI 
(right panel) shows the normalized MCE. Aside from the 
phase boundaries where again sharp peaks appear it is 
largely determined by the behavior of the saturation field 
(Fig. [3]) since the bare Tn^dh ^ h'^) is roughly constant 
in (f). 

Deep within the ordered phases the degree of enhance- 
ment of the MCE relative to an ideal paramagnet is 
chiefly controlled by the saturation field hc(4>)- The 
FTLM and LSW predictions are therefore in excellent 
qualitative agreement (c.f. Fig. [51 and Fig. [TH)) . However, 
once again, finite size effects prevent the FTLM method 
from capturing the full extent of the anomalous enhance- 
ment of the MCE in the highly frustrated regions at the 
borders of the CAE phase. These may in practice be 
overestimated by LSW theory, since it takes no account 
of new non-magnetic phases stabilized by fluctuations. 
None the less we can gain further insight into the strong 
enhancement of ri„c(^ ^ h^) and the singular peak at 
(j)c in Fig. 1161 within the LSW approach by expanding the 
spin wave energies of the fully polarized phase around the 
incipient ordering vector. Explicitly, for in the classical 
NAF sector we have 



= Ah + a2{ql + ql) - a^q^ + ql) + a'^q: 



(62) 



where q = k — Q is the distance from the NAF vector 
Q — (tt, tt) and Ah = h—hc is the excitation gap with the 
NAF critical field he = 8SJi, The expansion coefficients 
are given by 02 = 204 = S{Ji — 2J2) and a'^ — SJ2- 
This expansion may be inserted into Eqs. [^5]) and 
the integration performed approximately analytically. It 
is assumed that only modes with an energy Cq < fceT 
contribute appreciably to the integral. In performing the 
integration one has to distinguish two cases: If one is 
within the NAF sector the second order coefficient 02 is 
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FIG. 17: In this figure the subcritical fields /!.max(0) and ho{(j)) 
where Tmc{<i>) in Fig. 1141 (right panel) is maximal or changes 
its sign respectively are plotted as function of the frustration 
angle. In the CAP phase the characteristic fields are symmet- 
ric with respect to = O.Stt. 



nonzero. If one is at the boundary to the CAF regime 
02 = and the dispersion is determined by the mixed 
fourth oder coefficient a'^. 

Performing the integration in this (classical) limit 
k^T ^ Ah one obtains the approximate expressions 



po 



In 



kBT 



for 2J2 < Ji (NAF) 



he 

k^ 
for 2J2 



SJ2 



In 



/ kBT 



KiShAh) 
Ji (NAF/CAF) 



(63) 



(64) 



In the corresponding quantum limit k^T <^ A/j, both the 
heat capacity and the rate of change of magnetization 
with temperature have an activated behavior. However 
this cancels in between the numerator and denominator 
of Eq. lU to give 



h — hr 



(65) 



regardless of the degree of frustration present in the 
model. We note that this is exactly the form predicted 
at a quantum critical point on the basis of scaling argu- 
ments-4.. 

Returning to Eq. (|64l64p — in the first case, inside the 
NAF sector, the divergence in the MCE for h /i+ is 
only of a weak logarithmic type. However at the classical 
boundary to the CAF sector (which widens into the dis- 
ordered regime due to quantum fiuctuations) the singu- 
larity becomes a much stronger one essentially of inverse 
square root type. This is the reason that the MCE at 
h — h^ shows a large anomalous peak as function of (f> 
when crossing the NAF/CAF boundary. Because of the 
vanishing of the second order term in Eq. ([62]) the dis- 
persion has a saddle point at Q leading to a large DOS 
of low energy spin waves for 2 J2 ~ Ji and therefore a 
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stronger algebraic divergence of the MCE at he appears. 
The same arguments hold for the CAF/FM boundary. 

However a word of caution is appropriate. In our spin 
wave calculations we assumed that the classical mag- 
netic phases are stable throughout the phase diagram. 
Strictly speaking this is not true. As indicated in Fig.[l]in 
the shaded sectors around the CAF/NAF and CAF/FM 
phase boundary quantum fluctuations lead to instability 
of magnetic order and select a different nonmagnetic or- 
der parameter, presumably staggered dimer^^ and spin- 
nematioS.. While the broad features of our theory can be 
trusted, a truly quantitative theory for the nonmagnetic 
sectors would require to start from the proper order pa- 
rameter and their associated elementary excitations. In 
real materials, sufficiently close to hc{T), the critical be- 
havior of the MCE will also be sensitive to the details 
of interlayer coupling and magnetic anisotropy. Both re- 
finements remain as an outstanding challenge. 

It is also instructive to track the low to moderate field 
anomalies in rmc(/i, (/>) as function of (j) around the phase 
diagram. In Fig.[T7|the maximum field /imax('/') and the 
field ho{(j)) at which the MCE changes sign are plotted as 
function of frustration angle, normalized to the satura- 
tion field hc{(f>). While h^sx{4>) I hc{(t>) is rather constant 
throughout most of the range of angles, ho{(l))/hc{(f>) 
shows considerable variation in (j). The two characteristic 
fields are again symmetric around 0=0. Stt or Ji=0 for the 
same reasons as explained above. Note that the overall 
double-minimum structure of /io(0)/ hc(4>) as function of 
(j) prevents its use as a criteria to resolve the ambiguity 
of frustration angles that appears in zero-field thermody- 
namic considerations mentioned in Sec. IIIII It is obvious 
from the FTLM results in Fig. [10] (right panel), that the 
temperature dependence of the (normalized) MCE above 
he is strongly suppressed as h increases. This effect can 
also be understood from the LSW calculations. Approxi- 
mating Eqs. (|42l43p for small and large temperatures we 
obtain the ratio of the low and high temperature (nor- 
malized) MCE as function of h > he- 

fn,c(T« Jc) Ew^h^k) 

As long as the field is not too far above he there is still a 
considerable dispersion in e(/i, k) and the above ratio is 
larger than one (Fig.[Tni), i. e. T^c is T-dependent. Once 
h^ he however the dispersion is negligible compared to 
the gap energy A/i and the ratio in Eq. ([M]) approaches 
one, i. e. we recover the behavior of an ideal paramagnet. 

A comparison between FTLM and LSW predictions of 
the temperature dependence of the MCE for fields safely 
above he (in order to avoid the logarithmic singularity at 
he) is given in Fig. [THl There is a reasonable agreement 
in both magnitude and qualitative T-dependence. Note 
however that the LSW approximation becomes unreliable 
when T approaches Jc/fce and too many spin wave modes 
are thermally excited. 
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FIG. 18: Normalized MCE as function of temperature for 
above-critical field calculated with FTLM an spin wave 
(LSW) method. Frustration angles are <j!> = (NAF) and 
(j>/n — 0.74 (CAF) corresponding to the Sr compound. For 
T > Jc the FTLM result converges rapidly to 1. The LSW re- 
sult for finite h/hc has an asymptotic large temperature value 
different form one, only in the limit he it also approaches 
one. 



VII. SUMMARY AND CONCLUSION 

We have investigated the magnetocaloric properties of 
the J1-J2 model using the FTLM method for finite clus- 
ters and spin wave analysis starting from the classical 
magnetic structures. The one-magnon critical field or 
saturation field obtained from FTLM agrees well with 
the spin wave result. Finite size scaling results suggest 
that close to the CAF/FM boundary, the true critical 
field is determined by a two-magnon instability. This is 
consistent with the proposed existence of a spin nematic 
ground state in this parameter ranged. 

Both FTLM and spin wave results predict a strong en- 
hancement in the low temperature specific heat at the 
saturation field when the frustration angle crosses the 
phase boundaries. This may be explained by the large 
degeneracy of low lying states in these regions. At a con- 
stant intermediate field the specific heat exhibits a double 
peak structure around the NAF/CAF boundary. The en- 
tropy and specific heat show only moderate field depen- 
dence below the saturation field he. This feature may 
already be understood in a mean field approach where 
the entropy is strictly constant for all h < he- 

Likewise the strong enhancement of the MCE just 
above the saturation field was investigated. In the FTLM 
results, the MCE was enhanced by up to a factor ten rel- 
ative to an ideal paramagnetic (at temperatures small 
compared to the energy scale Jc). Surprisingly, the 
largest enhancement (from FTLM, relative to an ideal 
paramagnet) does not occur at the CAF/NAF boundary, 
but deep within the magnetically ordered sectors. This 
can be understood in terms of the anomalous enhance- 
ment of the specific heat in the frustrated regions, which 
enters into the denominator of the MCE (c.f. Eq. ([1])). 

The overall dependence of the MCE enhancement 
ratio reflects that of the saturation field. This is also 
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true for the spin wave results. There, in addition, the en- 
hancement is sharply peaked on the CAF/NAF boundary 
appears. This is due to the change of the field-scaling be- 
havior above the critical field from logarithmic to inverse 
square-root when the boundary is crossed. This feature is 
due to the appearance of Goldstone modes along a line in 
the BZ when has it critical value. The MCE enhance- 
ment may also be directly seen from the isentropics or 
adiabatic temperature curves which exhibit a large slope 
above the saturation field. 

Below the saturation field the results of the FTLM is 
strongly constrained by finite size effects, however the 
spin wave analysis provides considerable insight into the 
systematic properties of the MCE. In absolute terms the 
MCE is very small — in fact much smaller than the para- 
magnetic value. This is a direct consequence of mag- 
netic order, and can easily be understood from the mean 
field picture of a canted Neel state, for which entropy 
is constant as a function of field. The spin wave the- 
ory, however, predicts a flat maximum in the MCE at 
low fields, followed by a sign change for subcritical fields. 
These features are present throughout the phase diagram 
and the characteristic fields are moderately (^-dependent 
with a double-minimum structure. For this reason it is 
not likely that these low-field MCE features are useful in 
the determination of the frustration ratio. 

Considerably above the critical field the temperature 
dependence calculated from FTLM and spin wave theory 
show reasonable agreement. At temperatures of the order 
of Jc/fce the enhancement of the MCE is substantially 
reduced and the behavior of the system crosses over to 
that of an ideal paramagnet. 

The most pronounced discrepancies between FTLM 
and spin wave analysis appear at the classical phase 
boundaries. As explained above this may be well un- 
derstood in terms of the absence or presence of low ly- 
ing collective modes. While the former method under- 
estimates the MCE anomalies at the phase boundaries, 
the latter overestimates them — in fact it predicts a sin- 
gular behavior. A more advanced analytical treatment 
would have to take into account the quantum nature of 
the ground state around these boundaries, i. e. stacked 
dimer ( J2 > 0) or spin nematic ( J2 < 0) and the proper 
associated excitation spectrum. 

The present analysis provides some interesting predic- 
tions for the experimental investigation of the class of 
layered perovskites discussed in the introduction. Specif- 
ically we give detailed values for the possible saturation 
fields which should be easily accessible experimentally for 
Si04 and PO4 vanadates. According to Fig. 01 these crit- 
ical fields can be used to resolve the ambiguity in param- 
eterizing the model from its low-field susceptibility and 
heat capacity, far more cheaply than, e.g. neutron scat- 
tering. Furthermore we predict a genuine sign change in 
the MCE for subcritical fields which should be accessible 
to experiment. 

So far as practical applications — for example in cryo- 
gen free cooling — are concerned, the goal is to achieve 



as large a magnetocaloric effect as possible, at as low a 
field as possible. Here compounds not too far from the 
phase boundaries FM/NAF and CAF/FM are the most 
promising because they combine a significant MCE en- 
hancement with very moderate saturation fields. A de- 
tailed treatment of entropy as a function of magnetic field 
in the nematic phase occurring on the CAF/FM border 
remains an open challenge. However the high density of 
low energy excitations and low saturation field of this 
phase means that it looks a priori very promising for 
magnetothermal applications. 
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APPENDIX A: MOMENT DERIVATIVES 

Here we give the explicit expressions of the moment 
derivatives {d{Si) /d(3) {i =|| , _L) which are needed for the 
mean field calculation of the magnetocaloric coefhcient: 

{S^y = ^^ = [{I- A,,,)Bx+Ax,B,]/D (Al) 
op 

Where the determinant D is defined hy D = [(1 — 
Axx)i^~ ^zz) ~ ^xz-^zx] and the coefficients Ay {i,j ~ x 
or _L, z or II) and Bi are given by 



Bx = 
A — 

Axz — ^ 

and likewise 

B, = 



cosh^ i/?A 



(A3) 



— ^^(^1 tanhi/?A 

2A Acosh^i/JA A I A / 2^ 



2AAcosh^i^A A A2 2' 



jhz 



cosh^ i^A 



(A4) 



Azz — 



2A Acosh^ i/3A A \ A 



-^(l^] tanhi/3A 



"^^^ - 2AAcoshn/3A A A^ ^^nh ^/^A 



where a||,aj^ are given in Eqn. (|2ip and A,/i||,/ij^ are 
defined in Eq. 
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APPENDIX B: MATRIX FORM OF THE 
HAMILTONIAN 



bilinear in bosonic spin fluctuation operators. It can be 
written as 



In this appendix we give the expHcit matrix form of the 
Hamihonian in Eq. (I?!]) for the NAF structure which is 
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" A(/i,k) 
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(Bl) 



r 



where the constant term (arising from spin commutation 
relations) ensures that the zero point energy is negative. 
Using a simple coordinate rotation, we can reduce this 
matrix to a block diagonal form with two 2x2 diagonal 
blocks given by 



A{h,'k)±C{h,k) B{h,k) 

B{h,k) A{h,'k)±C{h,k) 



(B2) 



We can then solve each of the blocks using a separate, 
standard, Bogoliubov transformation to obtain the diag- 
onalized Hamiltonian in terms of NAF spin wave modes 
as used in Eq. ([55|l. 

Explicitly the complete transformation is given by 



axk = (ok + Aak) + I a[ 



A&LJ (B3) 



where A = ±1 is the branch index. The coefficients WAk 
and WAk of the transformation may be obtained by di- 



rect insertion and the requirement that the off-diagonal 
bilinear terms in the transformed operators a\k vanish 
identically. Alternatively the equations of motion for the 
ckAk may be set up and required to describe free motion 
with spin wave energy CAk- Both methods lead to the 
same condition on the coefficients given by 



2wAkWAk (^k + ACk) = (u^k + ^'Ak) ^k 



(B4) 



Using the representation UAk = cosh?7Ak, ^Ak = sinhr^Ak, 
one obtains the two branches (A = ±1) of the solution 
with 



1 



'7Ak 



■ tanh 



Bu 



Ak + ACk 



(B5) 



The prefactor of the remaining diagonal bilinear term in 
the transformed Hamiltonian gives the spin wave energies 
of Eq. 
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